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Abstract 

In this contribution we solve the following problem. Let H nm be a Hilbert space of dimension 
wn, and let A be a positive semidefinite self-adjoint linear operator on H nm . Under which condi- 
tions on the spectrum has A a positive partial transpose (is PPT) with respect to any partition 
Hn ® H m of the space H nm as a tensor product of an n- dimensional and an m-dimensional Hilbert 
space? We show that the necessary and sufficient conditions can be expressed as a set of linear 
matrix inequalities (LMIs) on the eigenvalues of A. 

1 Introduction 

This paper is motivated by a problem posed by Emanuel Knill and listed as Problem 15 in the list 
of open problems in Quantum Information Theory on the website of the Institute of Mathematical 
Physics at the TU Braunschweig jj. Knill posed the following question. 

Suppose A is a self-adjoint PSD operator on a Hilbert space H nm of dimension N = nm. Which are 
the conditions on the spectrum of A guaranteeing that A is separable with respect to any decomposition 
of H nrn as a tensor product H n Cg) H m of two Hilbert spaces of dimensions n, m, respectively? 

The problem was solved for the case n = m — 2 (a 2-qubit system) by Verstraete et al. [3], 
who showed that a necessary and sufficient condition is represented by a quadratic inequality on the 
eigenvalues of A. 

It is well-known that a necessary condition for separability is the so-called positive partial transpose 
(PPT) condition H|,[2|- Therefore any condition on the spectrum of A that is necessary for the PPT 
property to hold with respect to any decomposition H nm = H n (g) H m will also be necessary for 
separability with respect to any decomposition. 

In this contribution we explicite necessary and sufficient conditions on the spectrum under which 
A exhibits the PPT property (is PPT) with respect to arbitrary decompositions of fixed dimension. 
These conditions are expressed as linear matrix inequalities (LMIs) on the eigenvalues of A. Since 
the PPT property is equivalent to separability for the cases m — 2, n — 2, 3 01, we solve the original 
problem of E. Knill for these two cases and thus furnish an exact solution for a second special case. 
For arbitrary dimensions the obtained LMIs are necessary conditions for separability. The number of 
LMIs depends only on the minimum min(n, m) of the factor dimensions. 

Another result that emerged from our study is that the property of being PPT for arbitrary de- 
compositions of Hn into a tensor product of spaces of fixed dimensions n, m gets stronger if min(n, m) 
gets bigger. In other words, let N = nm = n'm' with min(n,m) > min(n',m'). If a PSD operator A 
is PPT with respect to arbitrary decompositions Hn = H n ® H m , then it is also PPT with respect to 
arbitrary decompositions Hn = H n i <g) H m i . 

The remainder of the paper is structured as follows. In the next section we derive the necessary 
and sufficient conditions on the spectrum of A in the form of LMIs. In Section 3 we show the above- 
mentioned dependence of the stringence of the condition on the dimensions of the decomposition. In 
the fourth section we illustrate our results for the cases min(n,m) = 2,3 and list the corresponding 
LMI conditions explicitely. In the last section we summarize our results. 

In the sequel H n denotes a Hilbert space of dimension n. Let A be a hermitian matrix of size 
nm x nm. We consider A as consisting ofraxm blocks of size n x n each. The partial transpose of 
A, denoted by A T , will be defined as an nm x nm- matrix consisting of the same blocks, but with the 
blocks (k, I), (I, k) interchanged for 1 < k < I < m. 
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2 Main theorem 



In this section we derive a set of LMIs on the spectrum of a self-adjoint PSD operator A on H nm which 
represents a necessary and sufficient condition for A to have the PPT property with respect to any 
decomposition of H nm as a tensor product H n £g) H m . 

Let x € R™ be an arbitrary real vector of dimension n. Let us define a set E{x) by 

E{x) = {x 2 k , (k = 1, . . . , n); +x k xi,-x k xi, (1 < k < I < n)}. 

This set contains n 2 real numbers. 

Our main result is based on the following lemma. 

Lemma 1 Let B — bb* be a hermitian PSD rank 1 matrix of size nm. Here b is a vector in H nm . 
Then there exists a real vector x € R p with non-negative, ordered entries x\ > xi > • • • > x p > 0, 
where p — min(n,m) 7 such that the spectrum of the partial transpose B T is given by the set E{x), the 
remaining p\n — m\ eigenvalues being zero. The elements x k , k — 1,. . . ,p are the singular values of 
the n x m matrix b that is obtained from b by arranging the m n-dimensional subvectors b%, . . . , b m of 
b columnwise. 

Remark: Obviously the assertion of the lemma holds also if B = 0. 

Proof. Let B — bb* satisfy the assumptions of the lemma, let b\ , . . . , b m be the n-dimensional 
subvectors of b and let b be the n x m matrix composed of these subvectors. 

Let now b = U n DV m be the singular value decomposition of the matrix b, where U n , V m are unitary 
operators of appropriate sizes and D a diagonal matrix of size n x m containing the singular values 
X\, . . . ,x p of b in decreasing order. 

Then W = U*®V m , Wv — U*®V m are unitary operators on H nm . It is easily seen that the relation 
(WBW*) T — WrB r W£ holds. In particular, the spectrum of B T equals the spectrum of (WBW*) r . 

Let us determine the structure of (WBW*) r . Let v\, . . . ,v m be the rows of V m . It is not hard to 
see that 

/ bx \ I U*hv* x \ 
w = Wb — (U* <x> V m ) : = : = vec(D), 

\b m ) \ U*hv* m ) 

where the operator vec stacks the columns of the matrix it is applied to into a big column vector. 
Hence the vector w — Wb has x k as (k + (k — l)n)-th element, k = 1, . . . ,p, and all other elements are 
zero. Partition the matrix WBW* = ww* into m x m blocks of size nxn. The block at position (fc, I) 
(k, I < p) has the product x k xi at position (fc, I), all other elements being zero. 

Let us consider the partial transpose (WBW*) r = (ww*) r of this matrix. Its block at position 
(k,l) has the product x k xi at position (l,k) (k,l < p), all other elements being zero. It is not hard 
to see that this matrix can be block-diagonalised by a permutation of the rows and columns, with 
diagonal blocks 

x2,(*=l,...,p); 7'),(l<*<Z<p), 

all other blocks being zero. Therefore the spectrum of the matrix (WBW*) T is given by x\ , k — 1, . . . ,p; 
+x k xi, —x k xi, 1 < k < I < p, the rest of the eigenvalues being zero. The first p 2 numbers form exactly 
the set E(x), where x is the vector composed of the singular values x%,...,x p of b. This completes 
the proof of the lemma. □ 

Lemma 2 Let n,m be positive integers and let p = min(n, m). Let further x 6 R p be an arbitrary 
vector. Then there exists a vector b £ H nm such that the spectrum of (bb*) r is given by the set E(x), 
the remaining p\n — m\ eigenvalues being zero. 

Proof. Define b elementwise as follows. Let x k be the (k+ (k — l)n)-th element of b, k = 1, . . . ,p, and 
let all other elements be zero. Let the matrix b, as in the previous lemma, be composed columnwise 
of the m n-dimensional subvectors 6i, . . . , b m of b. Then b is diagonal with x k as diagonal elements. 
Hence its singular values are the absolute values \x k \. Applying Lemma ^ we get that the eigenvalues 
of (bb*) T are given by the products \x k \ 2 , k = 1, . . . ,p; +\x k \\xi\, — \x k \\xi\, 1 < k < I < p, the rest of 
the eigenvalues being zero. But this is exactly the spectrum we claimed (bb*) r to possess. □ 

Lemma 3 Let A, B be hermitian matrices of size nxn. Let ax, . . . , a n and bx, ■ ■ ■ , b n be their eigen- 
values in decreasing order. Then 
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Proof. The scalar product (UAU*,B) is an analytic function f(U) of the unitary matrix U, 
which ranges over a compact set. Let Uq realize the minimum /* of this function. Then the first 
order extremality condition states that for any skew-hermitian matrix S we have {SUqAUq,B} + 
(UqAUqS*, B) = {S, [UqAUq , B]) = 0. In other words, the commutator [U AUq,B] is hermitian. Since 
the commutator of two hermitian matrices is always skew-hermitian, the matrices UqAUq and B must 
commute. But then we can diagonalize them simultaneously by conjugation with some unitary matrix 
V. Moreover, we have (U AU^ B) = (VU AU£V*,VBV*). 

Therefore the minimum /* is given by mm. ae s n Y^k=i a a-(k)bk, where a ranges over all permutations 
of the indices 1, . . . ,n. This minimum is attained at the inversion a* defined by c*(fc) = n + 1 — k. 
Indeed, let a £ S n and fc < I be such that cr(fe) < cr(l). Then a a a.-\ < CL a (l) an( I ^fc ^ ^i- It follows that 
(a CT ( fe ) - a CT (i))(6fe -h) > and hence a a ^ k )b k + a^bi > a a{l) b k + a a ( k) bi. Therefore an interchange 
of o~(k) and a(l) can only decrease the value of the sum X)fe=i a <r(fe)^fe- Performing this interchange 
consecutively for all pairs (k, I) for which k < I and cr(k) < a (I), we finally arrive at the inversion a*, 
regardless of the permutation a we started with. This completes the proof. □ 

Let A be a self-adjoint PSD operator on H nm . Let Ai, . . . , A nm be the eigenvalues of A in decreasing 
order and let p = min(n,m). Define also p + = p(p + l)/2, p_ = p(p — l)/2 and let S+ = {(k,l) | I < 
k < I < p}, S- = {(k, I) | 1 < k < I < p} be sets of index pairs with cardinalities accordingly. 
Note that we have E{x) — {x k xi \ (k, I) £ S + } U {— x k x\ \ (k, /) £ S_} for x E R p . In the sequel we will 
consider orderings of the sets S+, S-. We define an ordering of a finite set S as a bijective map a from 
S onto the set {1, 2, . . . , #5}, where #S is the cardinality of S. 

Let x £ R p be a vector with non-negative entries. 

Definition 1 We say that an ordering o~+ : S+ — * {!,... ,p+} of 5+ is compatible with x if for any 
two index pairs (k\, (k2, h) £ S+ such that a+(kx,li) < o + (fc2,/2) we have x kl xi 1 > x k2 xi 2 . 

We say that a pair of orderings (er + : 5+ — > {I, . . . ,p+}, cr_ : S- — f {1, . . . o/ the sets S+, 5_ 

is compatible with x if o~+ is compatible with x and for any two index pairs {k\, h), (k2, h) S S— such 
that a+(k%, 1%) < cr+(fc2, fa) (recall that 5_ C 5+j we have <j_(fci, Zi) < cr_(fc2, ^2)- 

Thus for a compatible ordering I is the image of the index pair (fc, I) £ 5± for which the product 
x k xi is largest, 2 the image of the pair for which x k xi is second-largest and so on. Note that for any x 
there exists at least one pair of orderings which is compatible with x. 

Corollary 1 The operator A is PPT for all decompositions of H nm as a tensor product space H n ®H m 
if and only if for any vector x e R p with non-negative and ordered entries x\ > X2 > ■ ■ ■ > x p > 
there exists a pair (a + , <r_ ) of orderings which is compatible with x such that 

^2 X nrn+l-a + (k,l)XkXl - ^ (k,l) x kXl > 0. (I) 

(k,i)es+ (k,i)eS- 

Proof. By definition, A is PPT for all decompositions of H nm if for all unitary matrices U the 
partial transpose (UAU*) T is PSD. This is the case if for all vectors b £ H nm we have 

b*(UAU*fb = ((UAU*) r ,bb*) = {UAU*, (bb*) r ) > 0. 

Let us fix b for the moment. By Lemma ^ there exist non-negative numbers x\ > ■ ■ ■ > x p such that 
the spectrum of (66*) r consists of the set E(x), the rest of the eigenvalues being zero. Note that, apart 
from the zeros, there are p+ non-negative eigenvalues x k xi, (fc, I) £ S+, and p_ non-positive eigenvalues 
—x k xi, (k, I) £ S-. If (cr_|_, c_) is any pair of orderings compatible with x, then we have by Lemma El 

mmb*(UAU*) r b = mm(UAU*, (bb*) r ) = ^ K m +i-a+(k,l)XkXi - K-(k,l) x kXi- 

(k,l)eS+ (k,i)eS- 

If the expression on the right-hand side is non-negative for any vector x with non-negative and ordered 
entries, then the expression on the left-hand side is non-negative for all vectors b. Hence in this case 
A is PPT for all decompositions of H nm . 

Let now x £ R p be a vector with non- negative entries, and (cr+,cr_) a pair of orderings which is 
compatible with x. Suppose that 

2_j ^nm+l-<j + (kJ)XkXl - ^ ^cr-(k,l)X k Xl < 0. 

(M)es+ (k,0es_ 
By Lemmas |21 and |3 there exists a vector b £ H nm such that 

mmb*(UAU*) r b = Km+i-a + (k,i)XkXi - ^ K-(k,i)XkXi < 0. 

(k,i)es+ (k,l)eS- 
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Let Uq be the unitary matrix that realizes this minimum. Then we get that UqAUq is not PPT. Thus 
if A is PPT for all decompositions of H nm , then inequality holds for any vector x with non- negative 
entries and any pair of orderings (a+, er_) which is compatible with x. But such pairs do exist for any 
vector x. This completes the proof. □ 

Let us transform expression Q). Let (er+, er_) be a pair of orderings. Define apxp-matrix A(<r + , cr_) 
elementwise as follows. 



A w (cr + ,<7_) 

Then we have 



\im+l-a + {k,l) i k<l, 
-K-(l,k)i k > I. 



Km+l-a+(k,l)XkXl - ^ K-(.k,l)XkXl = X T A(<7+ , <J-)x. (2) 

(k,i)es+ (k,i)eS- 

The matrix A has thus signed eigenvalues of A as elements, namely the p+ smallest eigenvalues in the 
upper triangular part, including the diagonal, and the p- largest eigenvalues with a minus sign in the 
lower triangular part. The particular arrangement depends on the orderings er + , cr_. Define the finite 
set of pairs 

E± = {(<r+,fT_) 3 x\ > X2 > ■ ■ ■ > x p > : (<7+,cr_) compatible with x — (xi, . . . ,x p ) T } . 

Lemma 4 Let x G R p be a vector with non-negative ordered entries X\ > ■ ■ ■ > x p > 0. Then there 
exists a pair (<t + ,ct_) 6 S± which is compatible with x. 

Proof. Let x satisfy the assumptions of the lemma. Then there exists a sequence of vectors x k 6 R p 
such that linifc-joo x k — x and the components of x k satisfy the inequalities x\ > x\ > ■ ■ ■ > x k > 
for all k. By the definition of £±, for any k there exists a pair of orderings (a k _ 1 a k _) G £± that is 
compatible with x k . Let (er+,er_) be an accumulation point of the sequence {(erf, <7_)}. Obviously 
cr_) is in E± and compatible with x. □ 

Now we are ready to prove our main theorem. 

Theorem 1 The operator A is PPT for all decompositions of H nm as a tensor product space H n ®H m 
if and only if for all (<r+, cr_) G E± we /iawe 

A(ct+,ct_) + A(ct + ,ct_) t h 0. 



Froo/. Suppose that A(cr + ,cr_) + A(cr + ,cr_) T ^ for all (cr + ,er_) G E ± . Then x T K(a + , cr_)x > 
for all a; G R p and for all (<7 + ,er_) G E±. From Corollary^] relation and Lemma^it follows that 
A is PPT for all decompositions of H nm . 

Let now (er + ,er_) G E± and x G R p be such that a; T A(a+, cr_)o; < 0. We shall show that there 
exists a unitary matrix U such that UAU* is not PPT. By Lemma [21 there exists a vector b G H nm 
such that the matrix (bb*) r has the spectrum E(x), the rest of the eigenvalues being zero. Denote the 
eigenvalues of (bb*) r by fix, ... , Hnm- 

The expression on the left-hand side of (J2J is of the form 

nm 
fe=l 

where cr is a permutation of the index set {1, . . . , nm}. There exist unitary matrices U, V such that 
UAU* = diag(X a(1) ,X a ( 2 ), K(nm)), V(bb*) T V* = diag(pi, . . . ,/j. nm ). Hence 

x t A(<t+, a-)x = (UAU*, V(bb*) r V*) = ((V*UAU*V) r ,bb*) < 0, 

and V*UAU*V is not PPT. This completes the proof. □ 

Thus we have expressed the sought necessary and sufficient conditions on the spectrum of A as a 
set of LMIs. 



3 Decompositions of different dimensions 

In this section we prove that if N = n\m\ — ri2TO2 with min(ni,mi) > min(ri2,m2), than any self- 
adjoint PSD operator A on Hn which is PPT for any decomposition — H ni ® H mi , is also PPT 
for any decomposition Hn = H n2 (3 H m2 . 

In the sequel we explicite the dependence on the dimension p of the sets S+, SL, E± defined in the 
previous section, i.e. we write S+(p), S—(p), E±(p). Our result is based on the following lemma. 
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Lemma 5 Let p > q and suppose that the elements of the vector y G R 9 satisfy the inequalities y\ > 
2/2 > ' • ' > Vq > 0. Le£ er_) &e a pair of orderings of the sets S+(q), S-(q) that is compatible with 
y. Then there exists a vector x G R p with elements satisfying the inequalities X\ > Xi > ■ ■ ■ > x v > 
and a pair of orderings (p+, pJ) of the sets S+(p), S—(p) which is compatible with x such that p + equals 
<7+ on the domain S+ (q) of definition of cr+ and p~ equals cr_ on the domain S- (q) of definition of 

<7_ . 

Proof. Let y and (er + ,er_) satisfy the assumptions of the lemma. Define the vector x as follows. 

2 

For k < q let Xk — yk- For k > q choose Xk such that ^ > x q+1 > x q+2 > ■ ■ ■ > x p > 0. 

2 

Since ^ < 1, we have x q — y q > — > x q+ \. Therefore the sequence {xk} is strictly decreasing. 
Further, let (k,l) G S+(q) and (k',l') £ S+(p) \ S+(q). Then we have V > q. It follows that Xk<x v < 

2 

xix q +i < Xi— = x q < XkXi. Therefore p+(k',l') > p+(k, I) for any ordering p + that is compatible 
with i. As a consequence, any such ordering maps S+(q) onto the set {1, . . . , #S+(q)}. Let p+ be an 
ordering of S+(p) which is compatible with x. Define another ordering p' + of S+(p) by 



,, n ( cr + (fc,Z), (k,l)eS+(q), 
P+[ ' >~ I P+(k,l), (k,l)eS + (p)\S + (q). 



Since a + is compatible with y, and x^ — yk for k < q, we have also that p' + is compatible with x. Now 
it rests to choose p'_ as the unique ordering of S_(p) that makes the pair (p' + , p'_) compatible with x. 
□ 

Now we are ready to prove the above-mentioned result. 

Theorem 2 Let N = riirni = nim-i with p\ = min(ni,TOi) > p 2 = min(ri2,m2). Let A be a self- 
adjoint PSD operator on Hpj which is PPT with respect to any decomposition Hjq = H ni (8) H mi . Then 
A is also PPT with respect to any decomposition Hjy = H n2 (g) H„ l2 . 

Proof. Assume the conditions of the theorem. Let (er + ,(7_) £ Yi±{p2). By the preceding lemma 
there exists a pair of orderings (p + , p_) G S±(pi) such that p+ equals ct+ on S+(p2) and p_ equals rr_ 
on S-{p2). It follows that the upper left pi x p2-subblock of the matrix A(p + ,p-) equals the matrix 
A(o+, a-). But A(p+, p-) + A(p+, p~) T y by Theorem^and the assumption on A. Hence we have 
also A(cr+, (7_) + A((T + , (T-) T >r 0. Application of Theorem ^ completes the proof. □ 



4 Examples 

In this section we illustrate the obtained results on the examples of 2 x n and 3 x n bipartite spaces. 

Let m = 2, n > 2. Then we have p = 2, p + = 3, p_ = 1. Let us find the set E±(2). If x G R 2 is 
a vector with elements xi > X2 > 0, then we have the relations x\ > x\X2 > x\ on the products x^xz, 
(k,l) G S+(2). Hence the only element of S±(2) is the pair (o+, er_) defined by 

f (1,1) ^ 1, 
cr+ : < (1,2) i ^ 2 , cr_ : (1,2) 1. 
[ (2,2) ^ 3, 

The corresponding matrix A is given by 

a / \ I ^2n A2n-1 \ 

Hence we obtain the following necessary and sufficient LMI condition: 

/ 2A2n A2n-1 — Al 

\ A2n-1 — Ai 2A2n-2 

Since Afc > by the positivity of A, and Ai > A2 n -i, this matrix inequality reduces to 

4A2nA2n-2 — (A2n-1 _ Al) 2 > <^> Al < A2n-1 + 2y/\2n^2n-2- 

We obtain the following corollary. 

Corollary 2 Let A be a self-adjoint PSD operator on the Hilbert space Llm- Let Ai, . . . , A2 n be the 

eigenvalues of A in decreasing order. Then A is PPT with respect to any decomposition of H2 n as a 
tensor product H2 <8> H n if and only if the inequality Ai < A2 n -i + 2y 'A2 n A2 n -2 holds. □ 
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For n — 2 and n — 3 this inequality takes the forms 

Ai < A 3 + 2y/X^, (3) 



Ai < A s + (4) 

It is well-known ^ that for the cases m = 2, n — 2, 3 the PPT condition is equivalent to separability. 
Therefore we obtain also the following corollary. 

Corollary 3 Let A be a self-adjoint PSD operator on the Hilbert space H± (He)- Let Ai,A2,A3,A4 
be the eigenvalues of A in decreasing order. Then A is separable with respect to any 
decomposition of H4 (Hq) as a tensor product Hi &) Hi (Hi ® H3) if and only if inequality (^) 
holds. □ 

Inequality is already well-known to be a necessary and sufficient condition for separability for 
any partition of H4 as Hi ®Hi and was presented as partial solution of Knill's Problem 15 in quantum 
information theory We present here inequality as partial solution for the case m — 2, n = 3. 

Let now m = 3, n > 3. Then p — 3, p+ = 6, p_ — 3. We shall now determine the set E±(3). If 
x\ > xi > X3 > 0, then we have 

x\ > x\Xi > m.ax.(x\, X1X3) > mm(xi, X1X3) > X1X3 > x\. 

However, we can have both x\ > X1X3 and ^1X3 > x\. Hence S±(3) consists of two elements, and the 
corresponding matrices A are given by 

(A 3n A3„_i A3„_3 \ 
— Al A3„_2 A3„_4 , A2 = 
— A2 —A3 A3„-5 J 

Thus we obtain the two LMIs 

2A3„ A3„_i — Ai A3„_3 — A2 \ / 2A3„ A3„_i — Ai A3„_2 — A2 

A3n-1 — Ai 2A3„_2 A3„_4 — A3 | >Z 0, j A3 n _i — Ai 2A3„_3 A3„_4 — A3 1^0. 

(5) 




A3n-3 — A2 A3„_4 — A3 2A3„_5 / \ A3 n _2 — A2 A3„_4 — A3 2A3„_5 



Corollary 4 Let A be a self-adjoint PSD operator on the Hilbert space i?3 n . Let X\, . . . , X^n be the 
eigenvalues of A in decreasing order. Then A is PPT with respect to any decomposition of H 3n as a 
tensor product H 3 ® H n if and only if linear matrix inequalities @) hold. □ 



5 Summary 

In this contribution we presented necessary and sufficient conditions on the spectrum of a self-adjoint 
positive semidefinite operator A on a Hilbert space H nm of dimension nm under which A has a positive 
partial transpose for any decomposition of H nm as a tensor product space H n ® H m , where H n ,H m 
are Hilbert spaces of dimensions n, m. These conditions have the form of linear matrix inequalities on 
the eigenvalues of A (Theorem^. We showed that if these conditions hold for a pair of dimensions 
(n, m), then they hold also for a pair (n',m') (where nm = n'm') whenever min(n',m') < min(n,m) 
(Theorem |2J. Hence the condition is most stringent if n « m. 

For the case min(n, m) = 2 wc reduced the LMI condition to a single inequality (Corollary For 
the cases m = 2, n = 2, 3 our conditions are necessary and sufficient also for separability of A with 
respect to an arbitrary decomposition of the underlying space H4 or Hq (Corollary|3I). While the result 
for n = m = 2 was known before, the result for m = 2, n = 3 is new and provides a solution of E. 
Knills open problem number 15 for another special case. 
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